The universal equation of state of a unitary fermionic gas 
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It is suggested that for a fermi gas at unitarity, the two-body bond plays a special role. We propose an equation 
of state using an ansatz relating the interaction part of the /-body cluster to its two-body counterpart. This allows 
a parameter-free comparison with the recently measured equation of state by the ENS group. The agreement 
between the two over a range of fugacity (z < 5 for a homogeneous gas, and z < 10 for the trapped gas) leads 
us to perform the calculations of more sensitive quantities measured recently by the MIT group. 



Feshbach resonance makes it possible to adjust the strength 
of the inter-atomic interaction in a neutral atomic gas. When 
the scattering length goes to ±oo, there is no length scale 
left other than the average inter-particle distance and the ther- 
mal wavelength (assuming a zero-range interaction). The 
gas is then termed "unitary" and its properties are univer- 
sal when expressed in appropriate dimensionless units at all 
scales whether the system is fermionic or bosonic |l]]. Re- 
cent accurate measurements by the ENS group 10, [51] and the 
Tokyo group J3l have confirmed the universal nature of the 
equation of state (EOS) of a gas of neutral fermionic atoms, 
and have given fresh impetus to its theoretical understand- 
ing S More recently, direct measurements by the MIT 
group [7] of the isothermal compressibility k, pressure P, 
and heat capacity Cy jNks for a unitary gas have revealed 
the superfluid transition at T C /T F = 0.167(13). 

In this paper, keeping in mind the fundamental nature of 
the two-particle bond at unitarity, we propose a description 
of the unitary gas as consisting of singlet pairs, in terms of 
which all higher order clusters are expressed. The resulting 
Equation of State (EOS) extends the agreement with the ENS 
data J2, [3] on the grand potential over a much larger range 
of fugacity z than expected. However, this description breaks 
down for z > 5 for the homogeneous gas (and z > 10 for 
the harmonically trapped gas). For the homogeneous gas, 
z = 5 corresponds to a temperature T/Tp = 0.22, below 
which the proposed EOS cannot be trusted. We calculate, 
with our higher virial coefficients, the pressure, compress- 
ibility and heat capacity of the homogeneous gas to compare 
with the MIT data 0]. The calculation of these quantities is a 
stringent test since they require higher moments of the virial 
expansion. We find that inclusion of the higher virial coef- 
ficients yields agreement with the MIT data for pressure and 
entropy down to T/Tp = 0.3, and the compressibility and 
heat capacity to T/T F = 0.6. 

To set the stage for the proposed universal EOS, we briefly 
recapitulate the virial expansion of a two-component interact- 
ing homogeneous fermi gas Jit]. The grand potential f2(/3, fi) 
is defined as fl = — rlnZ, where r = ksT = 1//3 and 
Z is the grand-canonical partition function. Furthermore 
fi = — PV, and may be expressed in a power series of fu- 
gacity z = exp(/3/i), where \i is chemical potential. The 
grand potential fi = — tZi(0) Ylb=i bi z > where Z\((3) is 



ter integral. For an untrapped gas in volume V, we have 
Z\(Ji) = 2(V/X 3 ), where spin degeneracy of 2 is included 
and A = (2nfi 2 (3/m) 1 / 2 is the thermal wave length. For a 
harmonic oscillator (HO) trap in three dimensions, Z\{J5) — 
2/(Tiu>/3) 3 . For a unitary gas, the cluster integrals 6;'s are also 
temperature independent in the high-temperature expansion. 
Subtracting from Q the ideal part of the grand potential 
we obtain the interaction part of the EOS as 



n - O (0) = -rZi(p) J2(Abi)z 



(1) 



1=2 



where A6; = 6/ — b[ . Note that b\ = b^'' = 1, and cancels 
out on taking the difference. 

Consider now the special role played by Ab 2 of the two- 
particle cluster at unitarity. In such a gas, the spin-up 
fermions have a tendency to pair up with the spin-down 
fermions because the short-range interaction potential is on 
the verge of producing zero-energy bound states. The Fesh- 
bach resonance being in the relative s-state, ensures the pair 
interaction to be operative only between singlet pairs. One 
finds that |@] Ab 2 = (2y/2) x §(AZ 2 ), where the factor 2^2 
arises from the CM motion, AZ2 is the relative two-body par- 
tition function, and the "suppression factor" | arises from the 
fact that only half of the N particles can interact in a spin- 
balanced two-component Fermi gas. Note that 1 1 Qfl at unitar- 
ity AZ2 = |, yielding A62 = -75 for such a system. What 
about the A6;'s for the /-body clusters that appear in Eq. ([TJ? 
Keeping in mind that the unitary gas may be looked upon as a 
system consisting of forming and dissolving two-body pairs, 
we conjecture that for the scale invariant system, the Abi for 
I > 2 should be expressible in terms of (A& 2 ) with an appro- 
priate suppression factor. Viewing a /-body cluster as one par- 
ticle interacting with the rest from a cluster of (I — 1) paired 
particles, we assume that the suppression factor is given by 
2^(1-1), where Af(i-i) = (Z — 1)(Z — 2)/2 is in general the 
number of pairs in a cluster with (1—1) fermions. Thus our 
basic ansatz is 



do) 



Abi 



(-)' 



(A62) 
2-^(1-1) 



I > 2. 



(2) 



For 1 = 2, M\ = 0, and Eq. (0 is an identity. The alternating 
sign (— )' in the above equation was put in to keep the num- 



the one-body partition function, and bi is the /-particle clus- ber fluctuation (AN) /N = J2i ' hz ^iz not to grow 
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to a very large value with z, where (AiV) 2 = iV 2 — A^ 2 is 
the number fluctuation, proportional to the isothermal com- 
pressibility jllh . A large value of the compressibility would 
lead to a vanishing monopole excitation which is a signature 
of instability IU2I1 . 

Although our description of the higher virial coefficients in 
terms of the second may seem to be very different from the 
conventional one, similar relationship between the third and 
second virial coefficients have been found in anyons which 
is also a scale invariant system lfl3l [f4ll . This is obtained 
by demanding that the divergences in the three-body clusters 
cancel by similar divergences in two-body clusters in the high 
temperature limit. A formal derivation for arbitrary I for the 
unitary gas appears to be non-trivial. 

With this ansatz, 



n-n^ = -TZ 1 (p)(Ab 2 )J2(- 



1=2 



2-^(1-1) 



(3) 



Experimentally |2,|5|,[l5|], it is the quantity h(Q = fi/fi (0) 
that is extracted, where £ = 1/z. This is given by 



h(0 = l + (Ab 2 )- 



nc°) 



(4) 



In a homogeneous gas with a spin-degeneracy of 2, 

2 
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Viln(l + ze~ l ) dt 



It is worth noting that using Eq. © with Ab 2 = l/y/2, 
we obtain A6 3 = -1/2-3, A6 4 = l/8\/2, Ab 5 = 
— 1/64\/2, etc. Numerically A63 is known t o gr eat accuracy, 
it was calculated up to 8 decimal figures in 1 1611 and has now 
been improved to 12 decimal figures lfl7ll . Our ansatz for the 
third virial coefficient differs from the numerically computed 
value in the third decimal and as such cannot be exact. How- 
ever, as we shall see the agreement with EOS data is unaf- 
fected by such fine differences in A63. It is also estimated J2] 
that A64 sa 0.096 ± 0.015, and is consistent with our predic- 
tion within the error bars. It should be mentioned that A64 as 
quoted in fnll is different sign and magnitude from J2l and 
our value. This however destroys the agreement with the data 
from the ENS group. 

Before confronting the experimental data, we note that 
for a gas trapped in a three-dimensional harmonic oscillator 
(HO), Eqs. © is modified to QH 



/ l (C) = l + (A6 2 )£ 



1=2 
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(5) 



and ft (0) = \ J °° t 2 In (1 + ze~*) dt. The additional sup- 
pression factor of l// 3 / 2 in Eq. (© was derived in llrjll as- 
suming a local fugacity in a HO potential. 

We are now ready to compare our predictions given by 
Eq. © for the homogeneous gas and Eq. © for the HO with 
experimental data. In Fig.Q] h(() — 1 for the homogeneous 
unpolarized gas, as given by our Eq. ©, is plotted against 
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FIG. 1. (Colour online) The function h(Q — 1 for the untapped 
unitary Fermi gas (from Eq. as a function of £. The crosses 
represent the experimental data presented by Nascimbene et al. J3l 
in their Fig. 3a. 



£. The crosses on the plot are the ENS experimental data as 
found by Nascimbene et al. J2l The authors quote that h and 
( are accurate to within 6 percent. It will be seen from this 
figure that the series given by our Eq. is in good agreement 
with the data down to £ rj 0.2. To put the agreement in per- 
spective the same data are plotted as a function of z = l/£ in 
Fig. [2] along with the behavior of the EOS including virial co- 
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FIG. 2. (Colour online) The universal function — 1 for the the 
untrapped fermi gas plotted as a function of the fugacity l/£. The 
2nd, 3rd, and 4th virial expansions are also shown. The 4th order 
expansion labelled Virial4p has A64 = 0.096. The experimental 
data are the same as in Fig.[T] 

efficients up to the fourth order, as was done by Hu et al. JH 
We see that such a truncated series could match the data to 
about z w 1.7. Our series (01 extends this to about z w 5. 
This also underlines the importance of higher-order virial co- 
efficients A6;'s for I > 4, despite their rapidly diminishing 
values. For the curve labelled Virial4p we set A64 = 0.096, 
the estimated value yfl, rather than 0.088 given by our ansatz. 
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Note that the h(() has been calculated in |6] within Pade ap- 
proximation and including up to A&3. Despite deviation from 
ENS data for £ < 1, they obtain surprisingly good agreement 
for energy and entropy per particle down to T/Tp = 0.16. 

We now turn to the ENS measurement for the trapped uni- 
tary unpolarized gas, as extracted by Hu et al. and com- 
pare with our Eq. (0. (See Fig. [3]) Here the convergence of 
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FIG. 3. (Colour online) Fhe universal function h(Q — 1 for fermions 
in an harmonic trap as a function of (, Eq. (O. The crosses repre- 
sented the experimental data presented in Fig. 6 of Hu et al. 

the virial series is faster as expected, and the agreement is re- 
markably good down to ( « 0.1, Figure [4] shows this clearly 
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FIG. 4. (Colour online) The universal function of trapped fermions 
— 1 as a function of fugacity. See Eq. {5}. Virial4p is the 4th- 
order virial expansion with A64 = 0.096. The experimental data 
are the same as in Fig. [5] 

when the truncated predictions from previous work are com- 
pared with our result. The range of applicability of the virial 
series (0 is now extended fourfold to z ps 10. 

It should be noted that the series given in Eqs. and (0 
converge for any value of z. However, the range of validity 
of the sum depends on the maximum value of I as seen from 
Fig. 12 and Fig.|4]for free gas and HO respectively. The first 



few terms make a significant difference but the importance 
of Abi for I > 8 is minimal even at z = 5. The series gets 
saturated by the first twenty terms which is denoted as "full 
expansion" in all the figures. 

Encouraged by the agreement with h(z) even deep into re- 
gion z > 1 where the normal virial expansion is not expected 
to work, we next compare our predictions for the recently 
measured data on compressibility, heat capacity, and pres- 
sure by the MIT group |7]. Following their notation we write 
(3P = fp(X)/X 3 , where P is the pressure, X = hiz = pf/, 
and fp(X) is the universal function given by 



/ 00 
fp(X) = 2 / 5/2 (exp(X)) + A6, cxp(; X) 



(6) 



1=2 



The first term in the bracket is the contribution due to the 



ideal fermi gas, and is the standard fermi-dirac integral 118 1. 
All thermodynamic quantities can now be expressed in terms 
of this universal function and its derivatives. Specifically we 
have pressure and compressibility normalized by their zero 
temperature values, Pq , kq, given by 



P = 



P_ = 5T fpjX) 
P a 2T F f p (X)- 



k 3T f p {xy 



(7) 



where T/T F = An /[Sir 2 f P (X)] 2/3 is the dimensionless tem- 
perature scale and the prime denotes a derivative with re- 
spect to X. The heat capacity at constant volume and en- 
tropy are given by C v /Nk B - ^M*l _ _ 



4 f' P (x) 



^f{p-l/R),S/Nk B = ||f©-ln(z). The above expres- 
sions allow one to calculate the relevant quantities either as a 
function of fugacity z or temperature T/Tp using the virial 
expansion given by Eq. © and compare with the respective 
experimental data of Ku et al. 

In the light of our earlier remarks (see Fig. 1), these com- 
parisons are limited to z values less than 4.95, corresponding 
to T/Tp > 0.22. Fig. 5 shows the variation of the pres- 
sure, entropy, and heat capacity as a function of T/Tp. The 
agreement with experimental data improves noticeably as the 
higher A&;'s are included. The agreement for the pressure 
and entropy hold to T/Tp = 0.3, indicating that the first mo- 
ments of the virial expansion are good. This is not the case 
for the second moments, however, as the plots for heat ca- 
pacity vs T/Tp shows. The theoretical plots start deviating 
appreciably from the data for T/Tp < 0.6. The same behav- 
ior is seen in Fig. 6 where the compressibility is plotted as 
a function of pressure (in reduced variables). It is interesting 
to note that despite these deviations, a peak in the compress- 
ibility of about the right magnitude appears in the theoretical 
curve, though at a higher value of P/Pq or T/Tp. Though 
tempting, we are reluctant to interpret this as indicative of the 
onset of superfluidity in view of the inaccuracy of the virial 
description in this range of temperature or pressure. 

We conclude that the high-temperature virial expansion, 
in conjunction with our ansatz, can match the EOS over a 
significantly larger range of fugacity, corresponding to about 
T/Tp « 0.3 for the homogeneous gas. Our ansatz (given by 
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FIG. 5. (Colour online) Reduced pressure (top), entropy (middle) 
and heat capacity (bottom) shown as a function of T/Tf for the 
untrapped unitary Fermi gas. The experimental data are taken from 
Ku etal. [J- 



Eq. (|2]i) resulted from the picture of a unitary fermi gas as a 
dynamic collection of singlet pairs, and assumed that (A62) 
determines the higher virial coefficients. The resulting suc- 
cess of this picture may point to some truth in this conjecture, 
and poses a challenge for deeper understanding. 
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FIG. 6. (Colour online) Reduced compressibility shown as a func- 
tion of reduced pressure for the untrapped unitary Fermi gas. The 
experimental data are taken from Ku et al. @]. 
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